Estimation of Semileptonic Decays of B c Meson to S-wave Charmonia with NRQCD 
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We study the semileptonic differential decay rates of B c meson to S-wave charmonia, r\ c and 
J/^, at the next-to- leading order accuracy in the framework of NRQCD. In the heavy quark limit, 
TTi;, — > oo, we obtain analytically the asymptotic expression for the ratio of NLO form factor to LO 
form factor. Numerical results show that the convergence of the ratio is perfect. At the maximum 
recoil region, we analyze the differential decay rates in detail with various input parameters and 
polarization for J/ip, which can now be checked in the LHCb experiment. 
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I. INTRODUCTION 



Hadron collider provides a large amount of data on 
B c events. Wherein the most easily identified decay 
modes to tag the B c are: fully reconstructed channel 
B c — > J/fyir and semileptonic decay channel B c — > 
J/^f£u(£ = e,n). The CDF Collaboration made the first 
observation of the B c meson by the semileptonic decay 
channel at the Tevatron fourteen years ago Latter, 
the DO Collaboration performed the same analysis in a 
sample of 210 pb _1 of the Run II data [2|. The cross sec- 
tion of B c production at the Large Hadron Collider(LHC) 
is larger than that at the Tevatron by roughly an order 
of magnitude, which reaches 49.8 nb at the center-of- 
mass energy y/s = 14 TeV [!, H[. This makes the ex- 
perimental study of the differential branching fraction 
of B c meson semileptonic decays to charmonium feasi- 
ble. We can also obtain the information of the Cabibo- 
Kobayashi-Mashkawa(CKM) matrix element in B c de- 
cays, especially V c b which is not well determined. 

Recently, the BABAR collaboration measured the 
partial branching fraction AB/Aq 2 in bins of the 
momentum-transfer squared, with 6 q 2 bins for £>° — > 
n-l + u and 3 q 2 bins for B° -> p~l + v @. They found 
that the partial branching fraction of B° — > tt~£ + v de- 
creases as q 2 increases, while for B° — > p~l + v process, 
the partial branching fraction increases first and then de- 
creases as q 2 increases. Actually, we know that all of 
the five form factors in above two decay channels at the 
maximum recoil region increases with q 2 , at the next- to- 
leading order(NLO) accuracy according to the light cone 
sum rules calculation @,0]- The decrease of B° — > tt~£ + v 
is caused by the phrase space, which counteracts the en- 
hancement from form factors. In this work, we try to 
make out whether this happens or not in B c semileptonic 
decays to charmonia. 
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There exist several approaches in calculating the B c 
meson semileptonic decays to charmonium. Some of 
them are: the light cone QCD sum rules [cj-lllj|. the 
relativistic quark model (13 . fl3j ] , the instantaneous non- 
relativistic approach to the Bethe-Salpeter equation [l4j . 
the non-relativistic constituent quark model[16|, the co- 
variant light front model [lj|, and the QCD potential 
model [Tf . 

Consider that the B c meson is constituted by two 
heavy quarks with different flavors, which masses are 
much larger than the Aqqd, analogous to the situa- 
tion of heavy quarkonium, the system turns out to be 
non-relativistic. Hence the relative velocity of heavy 
quarks within the B c meson is small, i.e. v -C 1, 
though bigger than the velocities of quarks in charmo- 
nium and bottomonium systems, and the non-relativistic 
QCD(NRQCD) formalism is applicable to the study of B c 
meson semileptonic decays to charmonia. In the NRQCD 
framework, the matrix elements of the concerned pro- 
cesses can be factorized as 



(J/^McT^^B.) ~ ^(0)B^(0)^ (nc) T" .(1) 

n=0 



Here, T M = 7 M (1 — 75), the nonperturbative parameters 
iP(0)b an d V'(0)j/^(t )c ) are the Schrodinger wave func- 
tions at the origin for be and cc systems, respectively. 
T" are hard scattering kernels which can be calculated 
perturbatively. 

The paper is organized as follows: In section II we 
present the definition for relevant form factors and work 
out the expressions of form factors in the NRQCD frame- 
work. The dependence of the NLO semileptonic differen- 
tial decay rates on q 2 is also obtained. In section III we 
study the theoretical uncertainty, and analyze the result 
in detail of the maximum recoil region. The last section 
is remained for conclusions. 
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II. FORM FACTORS AND SEMILEPTONIC 
DIFFERENTIAL DECAY RATES IN THE 
NRQCD FRAMEWORK 

The B c — > J/4>(r) c ) transition form factors, /+, fp, V, 
Ao, Ai, and A2 are normally defined as follows (l8j ] 

2 2 

(77 c (p)|c 7 ^|S c (P)> - / + ( g 2 )(P^+^- mB °~ m,7c g^) 



+/o(9 2 ) 



9 y 7 



(2) 



+ {m Bc + mj/^)Ai{q ){e 



Here we define the momentum transfer q = P — p. Note 
that /o(0) = /+(0) at the maximum recoil. 

It is straightforward to calculate those form factors at 
the tree level in the NRQCD. They read 



F LO (g2) _ lGV2C A C F 7r('Sz + l) a ^(0)B>(0) J/t j4) 



(q'-iz-l) 2 ) 2 ^) m 3 N c 



(J/^( V ,s*)\c^b\B c {P)) 



2iV(g 2 ) 
m Ba + mj/jf, ' 



(J/1>(p,e')\c>yi>>r 6 b\B c (P)) = 2mj^A Q {q 2 ) £ —^-q^ 



A 2 {q 2 



S-q 



m% c - m 2 ^ 



if) 
_l 
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(q 2 - {z - l) 2 ) 2 z 3 / 2 m 3 N c 
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16V2C a C f tt^ITT (4z 3 + 5z 2 + 6z- q 2 (2z + 1) + l) a^(0)fl>(0)j/t 
(q 2 - (z - l) 2 f z 3 / 2 (3z + l)m 3 b N c 



(6) 



LO 2 16V2C a C f ttVz~TT(3z + l)ij(0) B M0)j/^ 



{q 2 - (z - l) 2 ) z 3 / 2 m 3 N c 



(7) 



rL(l , 8V2C A C F ir^z~TT(-q 2 



3z 2 + 2z + 3) a s ^(0) Bc 4>(0) Vc 



(q 2 - (z- l) 2 y Z 3 / 2 m 3 N c 
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8V2C A C F TrVz~TT (9z 3 + 9z 2 + llz - q 2 {5z + 3) + 3) a,#) Bc #), c 
{q 2 - (z - I) 2 ) 2 z 3 / 2 {3z + l)m 3 N c 

I 



(9) 



where z = m c /mb. 

In Ref . [ljl , the form factors of B c transition to rj c 
with alternative parameterizations had been calculated 
at the NLO accuracy in the non-relativistic limit. There 
are three typical scales of the process, which possess the 
hierarchy of Aqqd <C m-c <C m-b- Taking into account 
the work of Ref. [2l[ , where the NLO form factors of B c 
transition to J/vI/ are calculated, we expand the ratios of 
the NLO form factors to the leading order(LO) form fac- 
tors at first order in z = m c /mb expansion in the heavy 
quark limit — > 00. And the asymptotic expressions of 
which are then obtained analytically, that can be found 
in the Appendix A. 

For light leptons(^ = e,fi), their masses mg can be 
readily neglected, hence the semileptonic differential de- 



cay rate of B, 

dT 
dq 



r\ c iv depending on q 2 reads 
1 2 it/, 12 



r\ c lv) 



G F \V cb \ 
192tt 3 t 



-A(<z 2 ) 3 / 2 [/ + (<z 2 )] 2 . (10) 



B,. 



Here, G F is the Fermi constant; V c b is the CKM matrix 
element; and X(q 2 ) — (m^ +wi 2 c — q 2 ) 2 - Am 2 B jn 2 a . For 
lepton r, its mass can not be ignored in the analysis, in 
which the form factors /n and /+ are both involved in. 
Thus /o can be measured via B c — > t] c tv t process while 
is obtained in B c — > rj c £v£ decay. 
By virtue of the NLO form factors, we can easily gain 
the distribution of NLO differential decay rate on mo- 
mentum transfer q 2 . To check the convergence behavior 
of the ratio of NLO differential decay rate to LO one, 
we select three sets of different values of z and scale q 2 , 
and illustrate parameter dependence in the first diagram 



TABLE I: Theoretical parameters for different sets, with 
renormalization scale [i — 4.8 GeV, the lifetime of the B c 
t(B c ) = 0.453 ps, and G F = 1.16637 x 1CT 5 GeV -2 [H, 
where mi,, m c and A are given in GeV, while 1^(0)1 are given 
in GeV 3/2 [2||,m. 





m b 


m c 


A 


IV>(o)k 


l^(0)|, o 


IV»(0)|j/« 


set 1 


4.8 


1.5 










set 2 


4.9 


1.4 


0.10 


0.3615 


0.283 


0.283 


set 3 


5.0 


1.3 











of Fig Q] The parameters in the three typical choinces 
can be found in Table HI wherein the Schrodinger wave 
functions at the origin for J/^Sf are determined through 
their leptonic decay widths at NLO. 

For B c — > r\JLv channel, at the maximum recoil point 
q 2 = 0, we obtain a value 4.67±g;|| x Kr 12 |V c6 | 2 GeV" 1 
for (TTOl) . which is larger than the value of 0.65 x 
lO- 12 \Vcb 
model [lj^ 

recoil point tend to disagree with what in Ref. [l6| . In 
NRQCD calculation, the form factors of B c to r\ c are obvi- 
ously enhancing with q 2 increase, than other approaches 
light cone sum rules and nonrelativistic quark model, and 
the trend is sharpening at NLO, which counteracts the 
decrescence from the factor of phrase space. 

For the channel of B c — > Jj^lv{l — e,/i), the decay 
rates in transverse and longitudinal polarization of vector 
meson J/ty can be formulated as 



2 GeV 1 obtained in nonrelativistic quark 
. Besides, the results away from the maximum 



dTz 

dq 2 



G 2 F X{q 2 f' 2 \V cb \ 



2 q 2 



1927r 3 m| 



(11) 



dTr G 2 F X{q 2 Y' 2 \V cb \ 2 q 2 
dq 2 



1927r 3 m| 



\H + (q 2 



\H-(q 2 



(12) 

respectively. Here the helicity amplitudes are expressed 
as follows: 

v (i ) t — x/ 9M/9 — A m ) 



H ± (q' 



Ho(q') 



A(g 2 ) 1 / 2 
1 



Xiq 2 ) 1 / 2 



(13) 



2m 



A(<? 2 ) 



m Bc + mj/q, 



A 2 {q 2 ) 



+ (m Bc + m J /s S ,)(m% c - m 2 ^ - q 2 )A 1 (q 2 ) . (14) 

While summing up the various polarizations, the semilep- 
tonic differential decay rate of B c — > J/^lv over q 2 is 
obtained 



dT (B 



G 2 F X(q 2 ) 1 / 2 \V cb \ 2 q 2 



J/^lv) = 

X (\H + (q 2 )\ 2 + \H_(q 2 )\ 2 + \H Q (q 2 )\ 2 ), (15) 



1927r 3 m| c 

,2\|2 



> 
o 
(3 



m/m =1.5/4.8 



m /m =1.4/4.9 



m/m =1.3/5.0 



> 
0) 
CD 




m /m =1.5/4.8 



m /m =1.4/4.9 



m/m =1.3/5.0 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 



FIG. 1: NLO differential decay rates for the B c — > ry c ^(up) 
and B c —¥ J /^lu (down), for different values of quark mass. 
The renormalization scale is chosen to be close to the bottom 
quark mass, i.e. [i — 4.8 GeV. In the figure, Asymp. NLO 
means expanding the ratio of NLO form factor to LO one 
at the first order in z — m c /m b expansion and in the heavy 
quark limit mt —5- oo. 



V/* 



— q 2 ) 2 — 4m B 1 



where X(q 2 ) = (m B 

Similar as B c — > r\JLv, the distribution of NLO dif- 
ferential decay rate on momentum transfer q 2 for B c — > 
J j^lv channel with three sets of different values of z 
is illustrated in the low diagram of Fig. [TJ At the 
maximum recoil point q 2 — 0, we obtain a value of 
3-21±g:|f x 10- 12 |K&| 2 GeV -1 for UT5J), which is larger 
than the value of 0.6 x 10 -12 |Vd,| 2 GeV -1 obtained in 
nonrelativistic quark model 16]. Except for the enhance- 
ment from the NLO K-factor and the NLO Schrodinger 
wave functions at the origin, the result in LO in NRQCD 
calculation is intrinsically bigger than what obtained in 
nonrelativistic quark model. 



III. 



THEORETICAL UNCERTAINTY 



In this section, the theoretical uncertainties of the NLO 
scmilcptonic decay rates in maximum recoil region are 
investigated in detail. It is found that the main uncer- 
tainties of the concerned processes have two sources, the 
heavy quark masses and the renormalization scale. In the 
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TABLE II: The NLO partial decay width for various g 2 . Herein, £ may be any kind of light leptons, the e or fi. We take 
m c /mb = 1.4/4.9, fi = 4.8 GeV in the evaluation. The uncertainty in the first column of the value comes from varying 
the ratio of quark masses m c /mj from 1.5/4.8 to 1.3/5.0; while the uncertainty in the second column comes from varying 
the renormalization scale /i from 3 GeV to 6 GeV. For J/*?, the partial decay widths in transverse(£*jj an d longitudinal 
polarization (en ) of J/9 are presented separately. 



bins of q 2 (GeV 2 ) 


< q 2 < 1 


1 < q 2 < 2 


2 < g 2 < 3 


3 < g 2 < 4 


4 < g 2 < 5 


Ar(£> c -> (W- 15 GeV) 


o nr -+l. 17+1.96 
O.UU_ 77_o .74 


n 7Q+1.78+2.37 
y - ' °-1.16— 0.89 


i 9 n +2.78+2.95 
lz -"-1.77-l.ll 


1 c 9+4.47+2.73 
1J - z -3.77-1.41 


90 n + 7 -57+5.06 
^"•"-4.36-1.87 


Ar(B c -» J/9(e\)lv) (HT 15 GeV) 


n 7n+0- 215+0. 159 
u - '"-0.141-0.061 


9 R/ |+0.95+0.60 
z - o4: -0.60-0.23 


r a/l +2.54+1.34 
51-0.51 


1 1 98+6.03+2.61 
H- z °-3.35-1.00 


On Q7+14. 14+4.90 
ZU.a 1 _ 7.16-1.87 


AT(B C ->■ J/9(s\)lv) (KT 15 GeV) 


fi HI +1.14+1.40 
"• ul -0.78-0.53 


7 Q7+1-93+1.84 
' -°' -1.27-0.70 


10 64 +3 - 38+2 - 51 

J.U.O^_ 2 72 - .95 


14 OR+618+3.56 
±1.»U_ 3 62 _i .35 


99 n7 +12.06+5.31 
ZZ.U ( _6. 45-2. 01 


Ar(B c -> J/Mv) (10~ 15 GeV) 


R 71 +1.35+1.56 
' -0.92—0.59 


10 C9+2.89+2.45 
lu -" z -1.88-0.93 


1 R /1Q+5.92+3.86 
1U.1»_ 4 .24-1.47 


9fi 94+12.21+6.18 
z "- zy: -6.98-2.35 


43.04± 2 6; 2 ™ 



TABLE III: The NLO partial decay width of processes B c — > r\ a rv T and B c 
recoil point is at m 2 . We take the value of m 2 = 1.776 2 GeV 2 . 



J/9tv t for various q 2 , where the maximum 



bins of q 2 (GeV 2 ) 


m 2 < g 2 < 4 


4 < g 2 < 5 


Ar(B c ->■ /; to t ) (KT 15 GeV) 


9 * Rn +0. 9245+0. 655 
z '™ u -0 .538-0 .241 


I 7 £9+8.42+4.70 

I I -" z _4.56-1.73 


Ar(B c -> J/9tv t ) (10~ 15 GeV) 


n S91 +0.375+0.194 
U.OZ1_ 213-0 .073 


R Q99+4. 017+1. 648 
U.»z^_ 2 . 107-0. 625 



evaluation, we vary the charm quark mass to c = 1.4 GeV 
by ±0.1 GeV, the bottom quark mass to& = 4.9 GeV by 
±0.1 GeV and the renormalization scale /j, = 4.8 GeV by 
i^GeV. 

For light leptons(£ = e, /i), we divide q 2 into five bins 
in maximum recoil region(0 < q 2 < 5 GeV 2 ) and cal- 
culate the semileptonic decay rates separately. The re- 
sults are presented in Table [TTJ We find that at small 
q 2 {0 < q 2 < 1 GeV 2 ), the longitudinally polarized J/* 
events dominate over the transversally polarized ones 
by a factor 8.5, and the difference will be reduced as 
q 2 increase. While for lepton r, we divide q 2 into two 
bins (to 2 . < g 2 < 4, 4 < g 2 < 5 GeV 2 ) in maximum recoil 
region. Here the physical mass of lepton r is taken to 
be to 2 = 1.776 2 GeV 2 [22j], and the results are shown in 
Table HTT1 



investigated seperatedly. The distribution in the maxi- 
mum recoil is found testable in the LHCb experiment, 
and in turn the NRQCD factorization will be also tes- 
tified. To be noted that in the minimum recoil region, 
the NRQCD factorization is spoiled by the infrared di- 
vergences and hence is not applicable to analysis of those 
processes. 
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Appendix A: The NLO B c to Charmonia transition 
form factors 



IV. CONCLUSIONS 

The NLO semileptonic differential decay rates of B c 
meson to charmomia are analyzed in detail with various 
choices of parameters. The uncertainties of patrial de- 
cay widths in different bins of momentum transfer q 2 are 
discussed. For B c — » J/^iu process, the partial decay 
widths for transverse and longitudinal polarizations are 



In this appendix, the QCD NLO B c to charmonium 
transition form factors are given at the first order in 
power of m c /mi,. For compactness of the expressions, we 



define 



m c /m, 



b, s 



and 7 



Besides, 



the form factors at maximum recoil point, i.e. q 2 = 0, are 
also presented, which are in agreement with what given 
in references [19, ]2u . 
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+24(1 - 2s) 2 (s(2s + 3) - l)Li 2 (l - 2s)(s - l) 2 + 6s(2s + 5) (2s 2 - 3s + l) 2 log 2 (s) 
+3(s(s(2s log(2)(2s(76s - 193) + 289) + 4s(-120s 2 + 369s + (8s 3 - 52s 2 + 90s - 43) log 2 (2) 
-406) - 84 log 2 (2) - 281og(2) + 747) + 92 log 2 (2) - 1101og(2) - 116) + 16 log 2 (2) 
+4(7 - 9 log(2)) log(2) - 3) + 3(s(2s(s(2s(4s(4 log(2)s - 6 log(2) + 6) - 28 log(2) - 69) 
+156 log(2) + 113) - 2(13 + 58 Iog(2))) + 72 log(2) - 5) - 8 log(2) + 1) log(s) 
-6 (2s 2 - 3s + 1) 2 (16 log(2)s - 22s - 2 log( 7 ) + 4 log(2) + 1) log( 7 ) 

-12 (2s 2 - 3s + l) 2 (2s 2 + s - 2) Li 2 (l - s))} , (A8) 
Ag LO (0) 1 a 8r l,„„ n , 2/i 2 , 10n/ „ /l, 2 , , 119 

+71og(2)log(z) - 2 1 log(z) + 7 log 2 (2) + l^i 2 !) 
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-C A ( - | log 2 (z) - log(2) log(z) - | log(z) 

9 log 2 (2) 31og(2K 
4 8 Z 1 ' 



24 
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